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Abstract: 


A benchmark experiment on Z$U slab samples was conducted using a deuterium—tritium 
neutron source at the China Institute of Atomic Energy. The leakage neutron spectra within energy 
levels of 0.8-16 MeV at 60° and 120° were measured using the time-of-flight method. The samples 
were prepared as rectangular slabs with a 30 cm square base and thicknesses of 3, 6, and 9 cm. The 
leakage neutron spectra were also calculated using the MCNP-4C program based on the latest 
evaluated files of PSU evaluated neutron data from CENDL-3.2, ENDF/B-VIII.0, JENDL-5.0, and 
JEFF-3.3. Based on the comparison, the deficiencies and improvements in $U evaluated nuclear 
data were analyzed. The results showed the following. (1) The calculated results for CENDL-3.2 
significantly overestimated the measurements in the energy interval of elastic scattering at 60° and 
120°. (2) The calculated results of CENDL-3.2 overestimated the measurements in the energy 
interval of inelastic scattering at 120°. (3) The calculated results for CENDL-3.2 significantly 
overestimated the measurements in the 3—8.5MeV energy interval at 60° and 120°. (4) The 
calculated results with JENDL-5.0 were generally consistent with the measurement results. 
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1. Introduction 


Nuclear data are widely used in almost all important fields of nuclear science and technology, 
such as nuclear power, nuclear medicine, and nuclear safety, and they are the link between basic 
research in nuclear physics and the application of nuclear technology''*!. The accuracy, integrity, 
and reliability of nuclear data influences the development of the nuclear industry!®”!, A benchmark 
experiment is an integral measurement of the parameters that can be presented with a few 
approximations and assumptions. This method is effective for verifying the accuracy of nuclear data 
[8-10] The objective of this approach is to determine the differences and deficiencies between the 
measurements and predictions of different evaluated neutron libraries. The discrepancies between 
the measurements and calculations can guide the evaluator to flexibly adjust the parameters that 
influence the calculated cross-section. 

Uranium is a radioactive nuclide that contains three natural isotopes (8U, 7°U, and 74U). The 
natural abundance of PSU is 99.27%, which is the highest among the three isotopes. Advanced 


nuclear energy systems, such as the fusion—fission hybrid reactor!!!'7!, accelerator-driven system!!?,"4] 
and fast reactor''>!, have the ability to breed nuclear fuel and transmute waste, and SU plays a 
crucial role in emitting relatively fast neutrons in these processes. SU has been one of the most 
significant components of the nuclear fuel cycle, including fuel enrichment facilities, power reactors, 
and spent fuel storage, and it constitutes more than 90% of the fuel in power reactors!'®!7!, The 
Collaborative International Evaluated Library Organization of the Organization for Economic 
Co-operation and Development Nuclear Energy Agency, an international organization that evaluates 
nuclear data, first focused on identifying discrepancies in the evaluation of a few high-priority 
isotopes, including **°U"®!, A series of measurements are being planned at the neutron time-of-flight 
(TOF) facility of the European Council for Nuclear Research and the Geel Electron Linear 
Accelerator TOF facility at the Institute for Reference Materials and Measurements!'*!.With rapid 
advancements in nuclear technology, more accurate nuclear data of ?8U are required. 

The international evaluated databases all contained neutron data for Z$U, and we investigated 
the evaluation of the CENDL-3.2 ?°!, ENDF/B-VIII.0 ?!!, JENDL-5.0 7!, and JEFF-3.3 I] libraries. 
Table 1 lists the evaluated information for the SU neutron data from the four libraries. The data 
type of 8U in the four libraries included cross-sections of various reactions, angular distributions, 
and double differential spectra. All other libraries, except CENDL-3.2, contained covariance 


information. 
Table 1 Evaluated data of n+?38U reaction 
Library CENDL-3.2 ENDF/B-VIII.0 JENDL-5.0 JEFF-3.3 
Evaluation time EVAL-MAY19 EVAL-DEC14 EVAL-NOV09 EVAL-MARI17 
G. Y. Tang, Z. M. Shi IAEA O.I to, N CENDAN DEN, 
. Y. Tang, Z. M. Shi, . Iwamoto, N. 
Author g i a EC-JRC Geel 
B. S. Yu, G. C. Chen Consortium Otuka, S. Chiba 
COLL. 
Energy range 20 MeV 30 MeV 200 MeV 30 MeV 
Calculation code FUNF03 EMPIRE-3.2.3 CCONE TALYS 


The FUNF code can analyze a set of neutron nuclear data with incident energies ranging from | keV to 20 MeV, 
including cross-sections of all reaction channels, angular-energy distributions of emitted neutrons, 
double-differential cross-sections of various emitted secondary particles, and gamma production: multiplicity, 
cross-section, angular distribution, energy spectrum, etc. 

Most neutron data of PSU in the CENDL-3.2 library were evaluated in 1990, when FMT was 
adopted as the main code, and the theoretical model of this code has not been reported publicly. 
Subsequent updates mainly included the following: (1) In 2004, the angular distribution of 
secondary neutrons and gamma production were replaced with calculations based on the FUNF03 
codel“, The contributions of direct inelastic scattering were replaced with the calculated results 
based on the ECIS95°! and DWUCK4"*! codes. (2) The gamma production spectrum was revised in 
2008. (3) In 2009, the energy and angular distributions of secondary neutrons from inelastic 
continuum scattering were replaced with results calculated using the FUNFO03 code. (4) In 2019, the 
delayed gamma data were evaluated and compiled again, and the resonance parameters were 
adopted from the ENDF/B-VII.1 data file. 


A few experiments have been conducted on 8U to identify discrepancies in the evaluation of 
238U with D-T neutrons. In 1979, the leakage neutron and gamma spectra from ?SU spherical shells 
with different radii were measured by Lusia F. Hansen at Lawrence Livermore National 
Laboratory””!; In 2010, Y. Nie performed leakage neutron spectra measurements on a SU sample (a 
10 cm square base with a 5 cm thickness) using D-T neutron source at the integral experimental 
facility in the China Institute of Atomic Energy (CIAE), at measured angles of 45° and 1350281, 
Subsequently, in 2016, new measurements were performed at 60° and 120° on 8U samples with 
geometric dimensions of 10 cm square area and 2, 5, and 11 cm thicknesses!?”!, 

In this study, three larger *U samples were used to improve the count rate of leaked neutrons 
and obtain more accurate experimental data. The experimental samples were shaped as rectangular 
slabs with a square base of 30 cm and thicknesses of 3, 6, and 9 cm, corresponding to 0.81, 1.62, and 
2.43 mpfs, respectively, for a 14.5 MeV neutron. During the TOF measurements, neutrons leaking 
from the rectangular slabs were measured using a BC-501A detector at measurement angles of 60° 
and 120° from energy intervals between 0.8 and 16 MeV, using pulsed D-T neutrons. The results 
were analyzed to examine the *°U data files of CENDL-3.2, ENDF/B-VIII.0, JENDL-5.0, and 
JEFF-3.3. A Monte Carlo code, MCNP-4C°!, was used to perform the measurements. 


2. Experimental Setup 


The measurements were conducted at a benchmark experimental facility built on the pulsed 
neutron generator of the CIAE. A series of measurements on the slab samples have been previously 
performed at this facility®'*4!. Fig.1 shows the layout of the experimental setup. 
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Fig.1 Schematic for measurements of leakage neutron at CIAE. 
2.1 Neutron source 
The neutron source was produced using a 400 kV Cockcroft—Walton accelerator at CIAE, 
which were generated at the tritium—titanium (T-Ti) target with a surface density of 1.2 mg/cm? via 
bombardment of a pulsed deuteron (D*) beam. During the experiment, the energy and average 
current were 300 keV and 20 uA, respectively, at a repetition rate of 1.5 MHz. The neutron pulse 
width was as narrow as approximately 2ns. The energy of the pulsed neutrons was approximately 


14.5 MeV, with a yield of 2x10° n/s at the target. An associated particle system at 135° was 
established to measure the yield of the neutrons, which majorly included a silicon surface-barrier 
detector (SSD). Although neutrons are produced through the T (d, n)*He reaction, alpha particles are 
produced at the reverse angle. The silicon detector was used to measure alpha particles in the AQ 


solid angle at angle gy to the incident D* beam, and ¢,(6,£,) is the neutron flux in the unit solid 


angle at angle @ to the incident D* beam. These parameters can be calculated using Eqs. (2) and 
(3)5536]; 


N 
f (9,84 )= 764 (1) 
mxr’ 
AQ=—; (2) 


where E4 is the average energy of deuteron ions participating in the reaction to produce neutrons in 
the T-Ti target, Na is the count of alpha particles measured at angle @, Aq is the corrective factor 
related to the energy of the D* beam and the emission angle of a particles. AQ is the solid angle of 
the silicon detector subtending the T-Ti target, r is the radius of the aperture to limit the solid angle 
in front of the silicon detector, and L is the distance between the aperture and the target. 

The relationship between the neutron flux in the unit solid angle at angles 6, ¢,(6,E,), 
expressed by Eq. (1), and the neutron yield is expressed as follows: 

GAO) 2018) 
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where Gr is the total cross-section of the D-T reaction, o(@) is the differential cross-section of the 

D-T reaction at angle 0. By substituting Eq. (1) into Eq. (3), the yield of neutrons can be obtained 

using Eq. (4): 

N, — N, x Áz x O ot (4) 
AQ xo (8) 
Here, the energy of the D* beam was approximately 300 keV, and Ea was approximately 147 

keV. The alpha particles were measured at 135° (@=135° ).When the neutron was emitted at 0° 

(0 =0° ), oo: was 3984 mb, o(0) was 336 mb/sr, and A, was assumed to be 1.263. The radius of the 

aperture was approximately 0.16 cm, and the distance between the limiting aperture and T-Ti target 


was approximately 90 cm. The expression, N, =N,, x1.508x10° , can be obtained. Therefore, as 


long as the number of ‘He particles is accurately measured, it is logical to calculate the number of 
neutrons. 
2.2 Samples 

The slab samples used in the experiment included one polyethylene sample and three SU 
samples. Polyethylene was used as a standard sample to test the reliability of the system; its purity 
was higher than 99.9%.The purity of the Z$U sample was higher than 99.37%. Further descriptions 
of the PSU samples are provided in Table 2. 

Table 2 Size, density, and composition of 8U samples 


Sample Size(cm) Density(g/cm*) Nuclide(mass ratio) 


ae 18.654 238U(99,3770), 235U(0.2015), C(0.159), N(0.019), 
2338Y 30x30x6 18.649 Mo+Zn+Pb+Bi+Ta+Cd+Ni+Fe+Co+W+B+Si+Mn+Cr 
30x30x9 18.868 +Mg+Sn+V+Ca+Cu+A1(0.2486) 


2.3 Collimator 

During the experiment, neutrons interacted not only with the $U sample, but also with other 
devices in the laboratory. These neutrons entered the detector after scattering and forming 
background counts. Shields and collimators were used in the experiment to minimize the influence 
of the background and scattered neutrons on the experimental results’. As shown in Fig.l, a 
precollimator consisting of iron, polyethylene, and lead was set between the experimental sample 
and concrete wall along the flight path of the neutrons. The other collimator was placed inside a 
2-m-thick concrete wall made of the same materials as the precollimator. A shadow bar was placed 
between the experimental sample and the precollimator along the flight path of the neutrons to 
eliminate the scattered background neutrons generated during the interaction between the source and 
the collimator. By using the shadow bar, precollimator, and wall collimator, the background 
spectrum (sample-out) measured during the experiment was significantly lower than the effect 
spectrum (sample-in), and the foreground—background ratio increased significantly (Fig.2). 
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Fig.2 TOF measurement results at 60° and 120° under sample-in and sample-out conditions. 
2.4 Electronics system and data processing 
A simplified diagram of the measurement circuit is shown in Fig.3. The leakage neutron spectra 
measurements were conducted using a scintillation detector comprising a BC-501A scintillator with 
a diameter and thickness of 5.08 and 2.54 cm, respectively. The BC-501Asignal was divided into 
two outputs. One output entered the analog-to-digital converter as the pulse height (PH), and the 
other was used to determine the TOF and pulse shape discrimination (PSD). The TOF was derived 
by considering the anode signal as the start-time signal, and the beam pick-up signal from the 
accelerator was used as the stop signal after an appropriate delay. The PH was used for energy 
calibration to determine the detector threshold. The PSD was used for n-y discrimination to 
eliminate y-ray events. 
Stilbene and BaF2 scintillation crystals were used to detect the TOF spectra of the source 


neutrons and TOF spectra of associated gamma rays, respectively. After adequate processing of 
these spectra, the time distribution of the pulsed neutrons was obtained. The SSD signal was the 


output of the silicon detector and was used to obtain the alpha-particle count using the associated 
particle system. Digital data were recorded by event from the CAMAC module and analyzed offline. 

In the experimental data processing, the data-taking software, Kmax, was used to obtain the 
initial TOF spectra of the leakage neutrons after selecting an appropriate PH threshold and neutron 
events via offline analysis. The initial spectra were then normalized to the neutron yields calculated 
based on the associated particle counts and BC-501A detector area. Hence, the final TOF spectra 
were in the neutron flux per unit neutron per unit area (cm). 
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Fig.3 Block diagram of measuring circuits. 


3. Monte Carlo simulation 


Simulations are an essential aspect of benchmark experiments. Problems existing in the 
evaluated data file can be identified by comparing the measurements and predictions of different 
evaluated neutron libraries. The simulated TOF spectra of the leakage neutrons were calculated 
using MCNP-4C, in which the modeling of a 14.5 MeV source and neutron detector is crucial for 
obtaining accurate results. 

3.1 Description of neutron source 

A description of the neutron source, which primarily includes the energy, angular, and 
pulsed-time distributions, is shown in Fig.4. One of the calculations was performed using the 
TARGET®! code, which was developed by PTB in Germany to express the energy and angular flux 
distributions of the neutron source. By inputting the values of the incident energy of the D” beam, 
T/Ti ratio, and T-Ti target thickness, which in this study were 300 keV, 1.8, and 1.2 mg/cm’, 
respectively, a calculated distribution that could be used in the MCNP simulation was obtained. The 
results are presented in Fig.4a. 
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Fig.4.Description of D-T source in Monte Carlo simulation: (a) angular dependence of energy spectrum; (b) 

pulsed-time distribution. 

The pulsed-time distribution is a critical input parameter used in the simulation of the TOF 
spectrum and is described by a TME card in the MCNP code. It was obtained through a fitting 
process of the neutron spectra from the stilbene scintillation crystal and the y spectra from the BaF2 
scintillation detector“!. The fitting results for the pulsed-time distribution are presented in Fig.4b. 
3.2 Description of neutron detector (BC-501A) 

The description of the BC-501A detector mainly includes the efficiency curve and size of the 
liquid scintillator detector. The efficiency of the BC-501A detector was determined by following 
these steps'*°!, (1) The relative efficiency of BC-501A was measured using a **Cf neutron source. (2) 
The absolute efficiency was calibrated using a D-D neutron, and the results were used to normalize 
the relative efficiencies. (3) The efficiency was calculated using the NEFF code"™"!, and the results 
were validated using the experimental measurement results presented subsequently. The light output 
function used in the calculation was calibrated with neutrons generated via the Be (d, n)!°B reaction 
using a 2x13 MeV Tandem accelerator at the CIAE. Fig.5 shows the shapes of the efficiency curves 
for different threshold values. The calculations are in good agreement with the measurements. 
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Fig.5 Detection efficiency curves of BC-501A liquid detector with different thresholds. 


Detector size refers to the area and thickness of the detector. The detector can be regarded as a 
surface detector because its thickness was only 2.54 cm, which is significantly shorter than the flight 
distance of 8 m. Therefore, a ring detector estimator with a diameter of 5.08 cm was used in the 
MCNP-A4C code. 

3.3 Simulation of MCNP program 

The benchmark model for the two measured angles used in the simulation is shown in 
Fig.6!*7!.When calculating the leakage neutron spectra transmitted at 60°, the material of the sample 
cell at 120° was described on a material card assigned to air, whereas the sample cell at 60° was 
assigned to air for 120° calculations. Both sample cells were assigned to air for background 
simulation. 

The MCNP program was adopted to simulate the experiment and test the latest data of >8U in 
CENDL-3.2, ENDF/B-VIUI.0, JENDL-5.0 and JEFF-3.3, with data files of other materials, such as 
targets, collimators, shields, and air taken from the ENDF/B-VUI.0 library. The ACE format data 
suitable for the MCNP code were processed using the NJOY99 programi]. The number of 
simulated histories was 10°, making the statistical uncertainty for the calculation slower than 1% for 
each time bin. Based on these comparisons, it was easy to determine the differences between the 
four evaluated libraries. 
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Fig.6 Geometry for MCNP simulation. 


E Polyethylene 


4. Results and Discussion 


4.1 Experimental system test on standard samples 

The leakage neutron spectra measurements on standard sample (polyethylene) with 
measurement angles of 45° and 60° were first conducted to verify the reliability of data acquired 
using this experimental system“, The polyethylene sample was a rectangular slab with a 30 cm 
square base and 6 cm thickness. 

The MCNP program was used to obtain the corresponding simulated results, and the data file 
was obtained from the ENDF/B-VIIL.0 library. The results were analyzed to examine the simulations 
after background deduction and return. The results are plotted in Fig.7, and the count values are 
listed in Table 3. The experimental results of the n-p elastic peak satisfactorily replicated the 
simulated results. This indicates that the processing of the experimental data is reasonable, and the 
system measurement data are reliable. 
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Fig.7 Comparison between calculation results for polyethylene at 45° and 60° and experimental results. 


Table 3 Results for standard sample (unit source neutron and unit detector area) 


Angle Measurement Simulation C/E 
45° (8.42+0.28)E-11 8.25E-11 0.981+0.032 
60° (1.03+0.03)E-10 9.93E-11 0.960+0.031 


4.2 Uncertainty analysis of experimental data 

Normalization of the initial experimental spectra can be performed not only based on the 
number of associated particles, but also based on the ratio of the area of the n-p elastic peak between 
the experimental and simulated results. The normalized coefficient, B!, can be determined using Eq. 
(5): 


e 


Ne NG (5) 


p 
where N* is the count of neutrons (n-p scattering) on polyethylene, Nọ is the count of neutrons 
(n-p scattering) on polyethylene calculated using MCNP (unit source neutron), and N? is the count 
of associated œ on polyethylene. 

The initial experimental spectra on PSU sample can be normalized using Eq. (6): 


No NE Ng x NxN 
N, B'xN? = NUxNe 


n 


(6) 


where NE is the count of neutrons on 7*U (unit time bin), and N7 is the count of associated @ 


on 738U. 

Most systematic uncertainties in PSU measurements can be eliminated by normalizing the 
experimental data with the standard sample, as expressed by Eq. (6), which includes uncertainties in 
the absolute efficiency of the BC-501A detector (3%) and the ambiguity of the solid angle for 
associated particle detection. Finally, the total uncertainty of the experimental results consisted of 
the remaining systematic and statistical uncertainties. The statistical uncertainty from the neutron 
count from 7*8U was better than 3% over 80% points. The others are listed in Table 4. 

Table 4 Uncertainties in measurements 


source Uncertainty (%) 


The count of neutrons (n-p scattering) on polyethylene 


i <0.5 
calculated by MCNP (unit source neutron) 
Statistical The count of neutrons (n-p scattering) on polyethylene <0.5 
The count of associated œ on polyethylene <0.5 
The count of associated a on $U <0.5 
. The relative efficiency of neutron detection <3 
Systematic ae . 
The ambiguity of scattering angle <l 


4.3 Results and analysis 

Leakage neutron spectra measurements were performed at 60° and 120° on slab °U using D-T 
neutrons. The measurements were compared with the predictions of the MCNP calculations using 
four evaluated neutron libraries. The comparison showed that differences existed between the 
simulation and measurement results. In particular, the simulation with the CENDL-3.2 library was 
significantly higher than the experimental results between 200 and 400 ns. 
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Fig.8 Comparison between calculated and measured leakage spectra on *8U sample with thicknesses of 3, 6, 

and 9 cm at 60° and 120°. 

The discrepancies between the calculated and experimental results are considered to be 
attributed to the reaction cross-sections associated with the production of secondary neutrons in 
different libraries. NDplot was used to plot the experimental and evaluated data in the neutron 
libraries. The ND-plot program is a simple and efficient software application for nuclear data!*>!. By 
inputting the proton number, mass number, incident particle type, energy, and other information 
about the target nucleus, experimental and evaluated data, such as the reaction cross-section, 
secondary particle angle distribution, secondary particle energy distribution, and secondary particle 
energy angle distribution, can be obtained. The secondary neutrons produced through the interaction 
between the D-T neutrons (approximately 14.5MeV) and”8U are presented in Fig.9, which generally 


consist of neutrons scattered via elastic and inelastic scattering, and neutrons produced through 
fission reactions (n, f), (n, 2n), and (n, 3n). 
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Fig.9 Energy spectra of secondary neutrons after interaction between14.5MeV neutronsand”*8U sample at (a) 60° 
and (b) 120°.The data was retrieved from the CENDL-3.2 library. 

Based on the different energy regions mainly contributed by different reactions, the neutron 
energy spectra can be divided into four energy regions: (1) 12 < En < 16 MeV, which consist of 
neutrons scattered by elastic (n, el) and inelastic discrete-level ((n, inl)d) reactions; (2) 8.5 < En< 12 
MeV, in which the neutrons are mainly from inelastic continuum scattering ((n, inl)C) and total 
fission reaction (n, f); (3) 3 < En< 8.5 MeV, which consists of emitted neutrons from (n, f) and 
production of two-neutron (n,2n) reactions; (4) 0.8 < En< 3 MeV, which consists of emitted neutrons 
from (n, f), (n,2n), and (n, 3n) reactions and multiple neutron scattering. 

The simulated and experimental results of the leakage neutron spectra were integrated 
according to the above four energy regions, and the C/E ratios of the different libraries were 
obtained by comparing the simulated and experimental results. The C/E ratios of the four nuclear 
data libraries are listed in Table 5. 


Table 5 C/E ratios of leakage neutrons integrated over four energy regions for four libraries 
Thickness 


Energy 


CENDL-3.2 ENDF/B-VIII.0 JENDL-5.0 JEFF-3.3 
/cm /MeV 
60° 
12-16 1.385+0.046 1.006+0.033 1.020+0.034  1.081+0.036 
8.5-12 0.937+40.034 0.855+0.031 1.111+40.040 0.929+0.034 
3cm 3-8.5  1.331+0.043 1.035+0.033 1.03340.033 1.060+0.034 
0.8-3  1.043+0.033 1.107+0.035 1.082+0.034  1.053+0.033 
0.8-16 1.123+0.036 1.080+0.034 1.06940.034  1.054+0.033 
12-16 1.378+0.044 1.030+0.033 1.030+40.033 1.094+0.035 
8.5-12 0.847+40.028 0.781+0.025 1.044+0.034 0.845+0.027 
6cm 3-8.5  1.328+40.042 1.035+0.033 1.003+0.032 1.011+0.032 
0.8-3  1.031+0.033 1.065+0.034 1.03940.033 0.995+0.031 
0.8-16  1.105+0.035 1.052+0.033 1.031+0.032 1.001+0.032 
12-16  1.297+0.043 1.032+0.034 1.000+40.033 1.051+0.035 
8.5-12 0.913+0.033 0.791+0.029 1.016+0.037 0.839+0.031 
9cm 3-8.5  1.313+0.042 1.015+0.032 0.973+0.031 0.981+0.031 
0.8-3  1.001+0.032 1.012+0.032 0.983+0.031 0.933+0.030 
0.8-16 1.076+0.034 1.010+0.032 0.982+0.031 0.947+0.030 
120° 
Som 12-16 1.909+0.089 1.203+0.056 1.026+40.048  1.384+0.063 
8.5-12  1.278+0.051 0.968+0.039 1.171+40.047 1.099+0.044 


3-8.5  1.331+0.042 1.038+0.033 1.0300.033 1.064+0.034 
0.8-3  1.003+0.032 1.066+0.034 1.04840.033 1.014+0.032 
0.8-16 1.082+0.034 1.061+40.034 1.0454£0.033  1.029+0.033 
12-16 1.828+0.076 1.19540.050 1.004+0.042  1.351+0.056 
8.5-12 1.326+0.049 1.018+0.038 1.165+£0.043  1.108+0.041 
6cm 3-8.5  1.355+0.043 1.064+0.034 1.049+0.033 1.054+0.034 
0.8-3  1.013+0.032 1.0570.033 1.038+0.033 0.999+0.032 
0.8-16 1.093+0.035 1.059+0.034 1.042+0.033 1.014+0.032 
12-16  1.8530.073 1.183+0.047 1.01740.040 1.295+0.051 
8.5-12  1.2750.047 1.043+0.038 1.149+40.042 1.067+0.039 
9cm 3-8.5  1.328+0.042 1.045+40.033 1.013+40.032  1.016+0.032 
0.8-3 = 1.010+0.032 1.042+0.033 1.023+40.032 0.977+0.031 
0.8-16 1.078+0.034 1.044+0.033 1.022+£0.032 0.989+0.031 


Based on the C/E ratios, some deviations were observed between the simulated and 
experimental results. 

(1) In the 12 < En< 16 MeV region, the neutrons from elastic scattering were predominant at 
60°, whereas at 120°, the neutrons from inelastic discrete-level scattering were predominant. The 
simulations with CENDL-3.2 yielded considerable overestimations in both directions, particularly in 
the direction of 120°, which reached approximately 80%. The simulations with JEFF-3.3 and 
ENDF/B-VIII.0 satisfactorily reproduced the measurements at 60°; however, in the direction of 120°, 
overestimations of approximately 20% and 30% were obtained, respectively. The simulations with 
JENDL-5.0 accurately predicted the experimental values at both angles. 

The angular distributions of secondary neutrons from (n, el) and (n, inl)d reactions after the 
interaction of the 14.5-MeV neutron with the?ŝU sample are shown in Fig.10a and 10b, respectively, 
and the cross-sections of different libraries at 60° and 120° are listed in Table 6.The (n, el) and (n, 
inl)d cross-sections at 14.5 MeV in the CENDL-3.2 file are larger than those from the other three 
libraries at both angles. In particular, the data for the (n, el) reaction reaches 7.86E-03 b at 120°, 
which is significantly higher than that in the JENDL-5.0 file (2.50E-03 b). The larger (n, inl)d 
cross-section (9.42E-03 b) in the ENDF/B-VIII.0 file results in the overestimation of the 
calculations at 120° in this region. The overestimation of the calculations with the JEFF-3.3 file at 
120° is principally because of the larger (n, el) cross-section (3.51E-03 b) in it. 
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Fig. 10 Angular distribution of secondary neutron from different reactions for™°U at 14.5 MeV, retrieved 
from four data files: (a) elastic scattering (n, el), (b) inelastic discrete-level scattering (n, inl)d, and (c)inelastic 
continuum scattering (n, inl)C. 
Table 6 (n, el) and (n,inl)d reaction cross-sections for PSU at 14.5 MeV obtained from four different data files 
(unit:b) 
library degree (n, el) (n,inl)d (n, el) + (n, inl)d 
60 9.35E-02 5.05E-02 1.44E-01 
120 7.86E-03 1.32E-02 2.10E-02 
60 5.92E-02 2.58E-02 8.50E-02 
120 2.24E-03 9.42E-03 1.17E-02 
60 5.85E-02 3.05E-02 8.90E-02 
120 2.50E-03 7.25E-03 9.75E-03 
60 6.52E-02 4.20E-02 1.07E-01 
120 3.51E-03 1.20E-02 1.55E-02 


CENDL-3.2 


ENDF/B-VIII.O 


JENDL-5.0 


JEFF-3.3 


(2) In the 8.5 < En< 12 MeV region, the simulations with the CENDL-3.2 file showed a good 
agreement with the measured spectra at 60° but yielded an overestimation of approximately 30% at 
120°. However, the simulations with ENDF /B-VIII.0 was consistent with the measurements at 120° 
but yielded an underestimation of approximately 20% at 60°. The calculations with the JENDL-5.0 
file somewhat replicated the measured spectra at 60° but gave a slight overestimation at 120°. The 
calculations with the JEFF-3.3 files lightly underestimated the measurements at 60° but accurately 
predicted the measurement results at 120°. 

The angular and energy distributions of secondary neutrons from the (n, inl)C reaction after the 
interaction of the 14.5-MeV neutron with SU are shown in Figs.10c and lla. The angular 
distribution and energy spectrum provided by the CENDL-3.2 file at 120° were lower than those of 
JENDL-5.0, but the calculations using CENDL-3.2 appeared to overestimate the production of 
neutrons more than the JENDL-5.0 calculations. This effect mainly resulted from the contribution 
by neutrons partly scattered by the (n, el) reaction in the 8.5-12 MeV energy region owing to the 
wider width of the pulsed beam from the measurements. The angular distribution of the secondary 
neutrons at 120° obtained by the CENDL-3.2 file was higher than those of the other libraries 
(Fig.10a). 
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Fig. 11 Energy spectra of secondary neutrons from several reaction channels for ?8U at 14.5 MeV, retrieved 
from four data files: (a) continuous-level elastic scattering (n ,inl)C, (b) fission reaction (n, f), and (c) (n, 2n) 
reaction. 

(3) In the 3.0 < En < 8.5 MeV region, the simulations with ENDF/B-VII.0, JENDL-5.0, and 
JEFF-3.3 libraries show good agreements with the measured leakage spectra. However, the results 
obtained with CENDL-3.2 significantly overestimated the measurement results at both 60° and 120°. 
The discrepancy between the measured spectrum and that calculated with the CENDL-3.2 file is 
mainly because the energy distributions of secondary neutrons from fission and the (n, 2n) reactions 
are harder than those in the other three libraries. As shown in Figs.11b and 11c, the shape of the 
energy spectrum of secondary neutrons from the fission reaction in CENDL-3.2was higher than 
those from other files in the 3-8 MeV energy interval but lower in the 0.8-3 MeV interval. 
Moreover, the spectrum from the (n, 2n) reaction in CENDL-3.2 was higher than those from other 
files in the 3.5—-8 MeV interval but lower in the 1-3.5 MeV region. 

(4) In the 0.8 < En< 3.0 MeV region, calculations performed using the four libraries accurately 


predicted the measurement results. 


(5) Across the entire energy region (0.8—-16 MeV), the calculation values obtained with the 
CENDL-3.2 file were slightly higher than the measurement values, and the C/E ratios were closer to 
1 with an increase in sample thickness. 

A comparison between the measured spectra and those calculated with the CENDL-3.2 and 
JENDL-5.0 files was performed at each time bin, and the C/E ratios were plotted (Fig. 12). The C/E 
ratios of the leakage spectra calculated with the CENDL-3.2 file to the measured spectra tended to 
deviate from 1 over the high-energy region; however, the ratios became closer to | with decreasing 
energy. The C/E ratios of the spectra calculated with the JENDL-5.0 file to the measured spectra are 
generally very close to 1 in most energy regions, except for those that are slightly higher than 1 at 
approximately 10 MeV. 
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Fig. 12 C/E ratios for leakage neutrons leaking from SU with transport calculations obtained 
usingCENDL-3.2 and JENDL-5.0 evaluated data files. 


5. Summary 


An integral benchmark experiment on a rectangular-slab SU sample irradiated from a D-T 
neutron source was conducted to test the evaluated nuclear data. The leakage neutron spectra 
measurements were conducted at 60° and 120° using the TOF technique. The experimental results of 
the neutron TOF spectra were analyzed via Monte Carlo calculations, which were performed using 
238U) nuclear data collected from CENDL-3.2, ENDF/B-VIII.0, JENDL-5.0, and JEFF-3.3 files. The 
measured and calculated spectra were compared based on the shape of the leakage spectra and the 
C/E ratios integrated over five different energy regions. The calculations with the JENDL-5.0 file 
satisfactorily predicted the experimental results; however, the calculations with CENDL-3.2 
overestimated or underestimated the measurements in different energy regions. These discrepancies 
were caused by the inaccurate angular or energy distribution of secondary neutrons in the evaluated 
data file. In summary, the angular distribution of secondary neutrons from the (n, el) reaction and the 


energy distribution of secondary neutrons from the (n, inl) and (n, fission) reactions in 8U in the 
CENDL-3.2 library should be improved in future studies. 
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